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Abstract 

An elliptic cylindrical wave expansion method by using Mathieu functions is developed to obtain 
the scattering field for a two-dimensional elliptic cylindrical invisibility cloak. The cloak material 
parameters are obtained from the spatial transformation approach. A near-ideal model of the 
invisibility cloak is set up to solve the boundary problem at the inner boundary in the cloak shell. 
The proposed design provides a more practical cloak geometry when compared to previous designs 
of elliptic cylindrical cloaks. 
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Recently, achieving the invisibility of objects has received increased attention. Designs 
based on spatial coordinate transformations have been proposed [HI2JE1SE1E1EIIE]- The 
inhomogeneous and anisotropic cloaks obtained from the spatial transformation bend the 
incoming waves around the cloaked region, so that the fields after emerging from the cloak 
are the same as if the incident waves had passed through the free space. The cloaking princi- 
ple was demonstrated experimentally in the microwave regime [J]. Cloaks designs have been 
reported and investigated mostly for spherical and circular cylindrical geometries primarily 
due to the simplicity of analysis for structures that possess radial and axial symmetries. 
Recently, elliptic cylindrical cloaks have been reported based on non-orthogonal coordinate 
systems, the cloak being described by permittivity and permeability tensors with non-zero 
off-diagonal terms and verification being done numerically by full-wave finite-element simu- 
lations [gun]. 

In this Letter, we will study the scattering for an elliptic cylindrical cloak by using a 
semianalytical approach that is based on Mathieu functions. By taking the advantage of the 
elliptic cylindrical geometry of the structure, we focus our analysis on the two-dimensional 
(2D) elliptic cylindrical cloak, because the wave equation can be simplified in comparison 
with the three-dimensional (3D) case. As in the circular cylindrical case [Hj, we introduce 
a small perturbation into the ideal cloak to avoid extreme values (zero or infinity) of the 
material parameters at the cloak's inner surface. 

First, let us look at the wave equation inside the elliptic cylindrical cloak. In terms of 
the rectangular coordinates (x,y,z), the elliptic cylindrical coordinates (u,v,z) are defined 
by the following relations [Til IT2] 

x = /cosh w cos t>, y = f sinhw smv, z = z, (1) 

with < u < oo, < v < 2tt, and / the semifocal length of the ellipse. The contours 
of constant u are confocal ellipses, and of constant v are confocal hyperbolas. The z axis 
coincides with the cylinder axis. Similarly to the circular cylindrical cloak [H|8j, a simple 
transformation (u', v', z') with 

u' = — ^ u + a, v' = v , z' = z (2) 

can compress space from the elliptic cylinder region < u < b into the elliptic shell a < 
u' < b, u' = a being the inner confocal ellipse. Note that u, u', a, and b are dimensionless 
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quantities. Let the position vector in the original system be written as r = xx + yy + zz 
and in the transformed system as f' = xx' + yy' + zz' . Following the notations in [T3], the 
scale factors Qj(j = 1, 2, 3) of the transformation [Eq. Q] are found to be 

_ | dr/du' | (cosh 2m — cos 2m') 1 / 2 

^ 1 = | df'/du' | ~ T(cosh 2m' - cos 2v') 1 ' 2 ' 

_ | dr/dv' | (cosh 2u — cos 2m') 1 / 2 

^ 2 = | df'/dv' | = (cosh 2m' -cos 2m') 1 / 2 ' 
I <9rV <9;z' I 

^3 = | dff ,/ dz , | = 1, where T = (b — a)/b. 

Finally, the permittivity and permeability tensor components for the cloak shell can be given 

as 

<V_/V_ T <v _ /y _ ^ e^_/V_ cosh 2m - cos 2v ' 
eo yUo eo Ho T' e yUo T(cosh 2m' — cos 2m') ' 

where air is assumed for the ambient environment and the interior regions. In the following 

we consider the transverse-electric (TE) polarized electromagnetic field (i.e., the electri- 
cal field only exists in the z' direction). One obtains the following general wave equation 
governing E z * field in the cloak's elliptic cylindrical coordinates 

d ( 1 dE ? ,\ d ( 1 dE* 
+ 



f 2 (cosh 2m' — cos 2m') I du' \/v du' J dv' \fi u > dv' 
1 d 2 E 7 , d 2 E 7 , 



\e a dz' 1 dt 2 J 

An exp(-iuit) time dependence is assumed, where uj is the circular frequency. If we substitute 
Eq. (H) for e z >, /V) and /v we find 

3 2 £ z , 1 <9 2 £ 2 , f 2 k t 2 , , 
T— — |- + — -4- + - — cosh 2m - cos 2m' = 0, 6 
du' 2 T dv' 2 2T v ; v ; 

where k t = ko — k z , ko = uj to/iQ, and fc 2 is the ^-component of the wave vector (for 2D 

case, k z = 0). Equation Q can be solved by a separation of variables E z > = S(v')R(u') and 

the introduction of a separation constant c, 

A 2 

dM' 

d 2 / , u'-a 

— o — c — 2q cosh 2 — — 
du' 2 V r 

where g = f 2 k 2 /A. Outside the cloak Eq. (|7| is the same, whereas Eq. (|8| is replaced by 

r d 2 i 

<^ — 2 - (c - 2g cosh 2m') ^(m') = 0. (9) 
l du' J 



r d 2 i 

\ — 2 + (c - 2gcos2M') = 0, (7) 
I dv' J 

T 2 ^ - (c-2gcosh2^— -) ^(m') = 0, (8) 



Equation |7| is known as the angular Mathieu equation [HI [121 E] • The solution is denoted 
by S pm (q, v, n) where p, m denote even (e) and odd (o), n denotes the order. For simplicity, 
in the following the dependence of q is skipped. Equation ^ is known as the radial Math- 
ieu equation. Similarly to the circular cylindrical coordinates where the radial solution is 
expressed in terms of Bessel functions J v (-), Y u (-), Hjp(-), and Hl 2 \-), the radial Math- 
ieu equation has four kinds of solutions [HI H2J [II] : J pm , Y pm , H pm i, and H pm 2, where 
p, m = e, o. Equation ^ transforms into Eq. ^ with the change of variable u' — *> In 
the following the primes of u' and v' coordinates are dropped for aesthetic reasons. 

Now, let consider a plane wave is incident on the cloak. In the 2D case, the incident field 
can be expanded in the u, v coordinates with the following expression [HJ [12] 

£*» = V^J2 inj pm(u,ri)S pm (v,n)S pm ((j),n)/N pm (n), (10) 

n 

where is the angle of incidence measured from the x axis in the (x, y) plane, and N pm is the 
normalization constant. Since Y pm has a singularity at u = a, in order to circumvent this, 
we introduce a small perturbation to the ideal cloak pj. We consider the inner boundary 
located at u = a + 5 where 5 is a very small positive number. Now the electric field in each 
region can be given by 



u 



> b) E z = a/87t i n J pm (u, n)S pm (v, n)S pm {<t>, n)/N pm (n) 



+ i n a p s ^(n)H pml (u, n)S pm (v, n)S pm ((/), n) 
(a + 5<u<b) E z = v^^Va^ {n)J pm {^,n)S pm {v,n)S pm {<f>,n) (11) 

n 

+ j^fiW^m i(^> n)S pm (v, n)S pm ((f), n) 
(u<a + 5) E z = \^^i n a { p l l (n)J pm (u,n)S pm (v,n)S pm (4>,n) 

n 

where a P m(n)(j = 1, 2, 3) are coefficients for the resulting field inside the cloak. The tangen- 
tial electric and magnetic fields E z and H v (which can be obtained from E z [IHE]) should 
be continuous across the interfaces at u = a + 5 and u — b, and the orthogonality of S pm 
allows waves in each order n and each combination p, m = e, o to decouple. One obtains 
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a W ( n ) = _ 



pm\ 

. (sc) (]= J pm (S/T,n) f P m{.a + ^n)-!3 pm {n)f pm (5/T,n) 
pm{ ] pm{ ] N pm (n)H pml (S/T,n)f3 pm (n)^ pm i(S/T jn )- J pm {a + 5,ny l "' 



a ( 3 ) / \ _ f3 pm (n)J pm (5/T,n) J^ pml (5/T,n) - ^ pm {bjT,n) 
pm N pm (n)J pm (a + 5,n) f3 pm (n)J^ pml (5/T,n) - / pm (a + i,n)' 

In these expressions we introduced the log- derivative function = F'/F, F being J pm or 
H pml , the prime denoting the derivative with respect to u, and 

_ C[S pm (v,n)]y[sinh\a + 6) + sin"v}^dv 
t[S pm {v,n)Y/[sm\i 2 {5/T) + sm 2 v]^ ( iv ) 

where the integration is performed numerically. Note that if 5 — > 0, then f3 pm (n) — > when 
= 0. Now, keeping almost the same notations like in [8], we obtain in the case of 2D 
circular cylindrical cloak for any order Z, 

(1) _ (in) 



a (2) = a (sc) = a (in) MkpS/T) JjMa + 5)} - pjijkoS/T) ( ;| 

;! ' °' //;'w/v) >'./// l; (/.>/V)-- /,/,•„(„ • ( sy' 



a ( 3) = Q ( m) PMhS/T) J% m (ko6/T) - JiihS/T) 
Ql ai J t [k Q {a + 5)} (3<%f\k 8/T) - ^[^(a + 5)] ' 

where aj is the coefficient of the incident plane wave, /3 = <5/[T(a + 5)], and the log- 
derivative function = F'/F now refers to the Bessel function J\ and Hankel function 
the prime denoting the derivative with respect to the argument. One can see the 



perfect similarity of the closed-form relations (12) and (14) for elliptic and circular 2D 
cylindrical cloaks. 

As an example, we consider a=0.8814, 6=1.4436, /=0.1m, and 5 = 10~ 5 . The semiaxes 
of the elliptic cylinder at u = b are 0^=0. 2236m, d y =0.2m, and at u = a, d x =0. 1414m, 
d y =0.1m. The wavelength of the incident plane wave in vacuum is A=0.15m. We found that 
n = 0-i-7 is sufficient for convergence of the calculated scattering field |14j . Figure [ija) shows 
the snapshot of the resulting electric- field distribution (i.e., the real part of E z normalized 
against the maximum absolute value) in the vecinity of the cloaked object. The scatterer 
is illuminated from the — x direction. The cloak's two confocal elliptic boundaries (inner 
and outer) are shown as black contours in the figure. The electric-field distribution clearly 



demonstrates the cloaking effect of the cloak shell to the incident plane wave. The maximum 
absolute value of the normalized field inside the cloaked object is 7.5 x 1CT 3 . There is still a 
little bit of the field perturbation on the back side of the cloaking shell. This perturbation is 
diminished at smaller values of the semifocal distance /, that is, when the elliptic cylinder 
approaches a circular one. An example is shown in Fig. [lib), where a=2.9982, 6=3.6895, 
and /=0.01m. Now, the semiaxes of the elliptic cylinder at u = b are d x =0. 2002m, d y =0.2m, 
and at u — a, d x =0. 1005m, d y =0.1m. The maximum absolute value of the normalized field 
inside the cloaked object is 1.5 x 10~ 4 . 

In conclusion, a semianalytical wave expansion method using Mathieu functions was ap- 
plied to study the electromagnetic scattering properties of a 2D elliptic cylindrical invisibility 
cloak. A small perturbation into the ideal cloak was set up to solve the boundary problem 
at the inner surface of the cloak shell. Closed-form relations for the expansion coefficients 
of the resulting field inside the cloak were given which are perfect similar to those obtained 
for the circular cylindrical cloaks. Besides the theoretical interest for the variable separation 
problem solved by using Mathieu functions in case of elliptic cylindrical invisibility cloak 
shells, the results could be also useful for designing this type of invisibility cloak due to the 
very simple permittivity and permeability tensors, in spite of the non-perfect invisibility of 
the cloaks having elongated elliptic shapes. 
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FIG. 1: Snapshots of the resulting electric field distribution in the vecinity of the cloaked object. 
The semiaxes on (x, y) directions of the outer and inner ellipses are (0.2236m, 0.2m) and (0.1414m, 
0.1m), respectively in (a), and (0.2002m, 0.2m) and (0.1005m, 0.1m) in (b). The incident light 
wavelength is 0.15m. 



